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In this proceeding, I give an overview of the latest progress in lattice-QCD calculations of the 
Roper resonance, and Roper-nucleon electromagnetic and axial transition form factors using dy- 
namical gauge ensembles. 
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Introduction: The nature of the Roper resonance, iV(1440) Pu, has been the subject of interest since its discovery 
in the 1960s. The mass of this positive-parity excited state is lower than the negative-parity ground state iV(1535) Sn, 
which is quite unusual, never being seen in mesons. This has led phenomenologists to postulate that the Roper might 
not simply be a radial excitation of the nucleon. Rather, it might be a hybrid state that couples predominantly to 
currents with some gluonic contribution or a five-quark (meson-baryon) state. Experiments at Jefferson Laboratory, 
MIT-Bates, LEGS, Mainz, Bonn, GRAAL, and Spring-8 offer new opportunities to understand in detail how nucleon 
resonance properties emerge from nonperturbative aspects of quantum chromodynamics (QCD). 

On the theoretical side, many attempts try not only to reproduce the reversal of the mass ordering of Roper and 
Sn starting from basic QCD properties, but also to determine the Roper-nucleon transition form factors. Among 
these theoretical efforts, lattice QCD (LQCD) is the most promising candidate for studying the nonperturbative 
aspects of QCD, where the Roper mystery lies. Lattice QCD has been successfully employed to study various strong- 
interaction phenomena, such as the spectroscopy of heavy-quark hadrons, the tower of excited baryon states, flavor 
physics involving the CKM matrix, hadron decay constants and the strong coupling constant. In many cases, it can 
provide higher-precision data from the Standard Model than what can be measured experimentally. By discretizing 
spacetime into a four-dimensional grid and working in Euclidean space, we are able to compute the path integral 
(in terms of discretized versions of the QCD Lagrangian and operators) directly via numerical integration, providing 
first-principles calculations of the consequences of QCD, without making uncontrollable approximations. The four- 
dimensional lattice size, lattice spacing, the quark masses in the sea sector and precision of the integration have been 
greatly improved over the past decade due to increasing computational resources and more efficient algorithms. We 
can systematically study the effects due to having a finite volume and nonzero lattice-spacing and remove them by 
extrapolation, using calculations with differing values of these quantities. 

Due to limited computational resources, many existing calculations have had to use multiple unphysical pion masses 
and extrapolate to the physical pion mass via chiral effective theory. (This helps to determine the low-energy constants 
in the theory.) The uncertainties caused by the extrapolation have been greatly reduced as lighter pion masses are 
included in calculations. Recently, major collaborations, such as BMW, PACS-CS and MILC have begun generating 
QCD gauge ensembles at the physical pion mass, which will remove such systematics and put stringent tests on low- 
energy effective theories. Lattice calculations with sub-percent precision have been achieved for many flavor-related 
quantities and used to guide tests of the Standard Model and to probe possible new physics. Calculations of hadron 
structure and interactions require more computational resources but should soon enter an era of precision. 

Resonance: Much progress has been made on the Roper resonance. The left-hand side of Fig. [T] was a state- 
of-the-art summary of lattice calculations of the Roper resonance as of the last Nstar workshop. Many groups have 
developed their analyses using the cheapest available lattices (having sufficiently large volume, fine lattice spacing, 
etc). However, most of these calculations have been done in the "quenched" approximation, where the effective sea- 
quark mass is infinite; that is, only gluonic interactions with the valence quarks are taken into account. Such an 
approximation preserves certain QCD properties, such as confinement, asymptotic freedom and spontaneous chiral 
symmetry breaking, making quenched simulations useful for testing new lattice techniques. However, quenching 
introduces uncontrollable sources of systematic error that are expected to be 30% or larger, depending on the physical 
quantity in question. 

In recent years, newer calculations have not needed to apply this approximation, and more results with Nf = 2 + 1 
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FIG. 1: (left) Summary of published Nf — LQCD calculations of the nucleon and Roper masses in GeV as functions of M%. 
Note that the errorbars are only statistical, (right) Summary of the latest Nf =2 + 1 Roper results as a function of pion mass 
squared. Note that the filled symbols and solid errorbars (open symbols and dashed errorbars) indicate the results taken from 
published papers (the latest Lattice Conference proceedings). 



(degenerate up and down plus strange in the sea) have become more widely available. 

Two preliminary studies were reported in the annual Lattice Conference using Nt = 2 + 1 domain- wall fermion and 
clover lattices by RBC/UKQCD and xQCD; their results are shown as the open symbols in Fig. [I] but no refined 
numbers have been reported since. The CSSM grouppQ used local baryon operators with different smearing parameters 
to construct a large basis for the nucleon with spin 1/2 and analyzed the data using the variational method. They 
use Nf = 2 + 1 isotropic lattices at a fixed lattice spacing, a = 0.0907 fm and fixed volume, 32 3 x 64 with lightest pion 
mass 156 McV. Their calculation of the Roper mass is shown on the right-hand side of Fig. [I] as 5 purple diamond 
points. However, the lightest two pion masses, since their M^L values are less than 3, may suffer finite-volume effects 
(which have not currently been estimated). Hadron Spectrum Collaboration (HSC) reported two studies using a large 
basis of independent operators organized according to representations of the cubic group, including operators where 
the three quarks inside the nucleon are located at different lattice sites[2j[3]. They have used 2 + 1-flavor anisotropic 
clover lattices, with finer temporal lattice spacing to improve the excited-state signal in Euclidean time; that is, the 
lattices have better resolution in the direction from which the particle energies are extracted. The green pentagons in 
right-hand side of Fig. [I] are taken from the lowest two states reported in their study last year using only operators in 
the Gig irrep. This year, a refined study mapped the cubic-irrep operators to definite spins using derivative operators 
constructed with continuum symmetries. They found four states closely aligned and roughly around the excited state 
found by CSSM. A refined CSSM study reported in this year's Nstar Workshop 1 using a larger basis of smeared 
operators that they found a similar structure as Ref. [3|. This may suggest that multiple-particle operators need to 
be included in the correlator matrix to clearly identify the nature of these closely related states; both HSC and CSSM 
are currently investigating this. 

The dynamical results in Fig. [I] are impressive, but the errorbars shown above are only statistical; that is, there 
is no estimation of the systematics due to the finite volume and nonzero lattice spacing in the calculations. This is 
sometimes difficult to do due to the limitations of the dynamical lattices available. Ref. |U has demonstrated with 
high-statistics calculations using Nf = 2 + 1 anisotropic clover lattices at M n ss 390 MeV. Four volumes are used in 
this calculation: 2, 2.5, 3 and 4 fm with M^L « 3.9, 4.8, 5.8, 7.7. The volume dependence is not noticeable in the 
earlier HSC study[5] with less trajectories available. To our surprise, the nucleon ground state with M n L around 4, 
which is usually the rule of thumb for where finite- volume effects can be ignored, shows around a 30-MeV discrepancy 
from its value in the infinite- volume limit. 

We investigate the effects of lattice volume on the Roper mass and found the effects to be even more dramatic; we 
observed the central- value shift to be about 170 MeV. More time should be devoted to studying volume dependence to 
get a better estimation of the finite- volume effects in these dynamical calculations. Furthermore, we found the overlap 
factors associated with the Roper states at the largest two volumes, gave a ratio Ar(A im)/An(3 fm) = 0.964(52), 
consistent with one; this suggests that the first-excited state is likely to be a one-particle state rather than a scattering 
state for the ensembles with pion mass larger than 390 MeV. 

Electromagnetic Transition Form Factors: Roper-A^ transition form factors provide a different approach to 
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FIG. 2: Proton-Roper transition form factor F[ R (Q 2 ) (left) and F% R (Q 2 ) (right) on the Nf = 2 + 1 anisotropic lattices with 
M w ~ 390, 450, 875 MeV whose volumes are 3, 2.5, 2.5 fm, respectively. 



understanding the unusual nature of the Roper state. They also provide information about the spatial distribution of 
the hadron. The general spin-1/2 electromagnetic form factors can be obtained from the following matrix elements: 

<JV a JViXqt) = u N2 (p') 

where q = p' — p. On the lattice, we can obtain both the form factors related to the Roper decay P\\ — > jN, and the 
one that is related to photoproduction, 7*7V — > Pn. 

We performed the first LQCD calculation of the Roper-nucleon transition form factors a few years ago, also starting 
with the quenched approximation where the QCD vacuum only contains gluons, using single pion mass around 
720 MeV; see Figs. 1, 2 in Ref. [§]. Since the nucleon and Roper masses are much higher than the physical ones, this 
calculation is in the time-like region for the matrix element (Pu\V^\N). When we decrease the pion mass, we will 
enter the space-like region and this matrix element will be helpful in giving us different Q 2 points. We found there 
are significant differences from the experimental values in the low-Q 2 region, which might be caused by the heavier 
pion mass used in the calculation or by the quenched approximation. 

A later study explored the pion- mass dependence of the transition form factors F* 2 at three pion masses, as well 
as large-Q 2 form factors; see Fig. 4 in Ref. [7J. In the time-like Q 2 regions, all four transition form factors have very 
mild quark-mass dependence. Even in the space- like region, the dependence only becomes evident at Q 2 > 2 GeV 2 . 
We obtain form factors with momentum transfer up to 6 GeV 2 ; in the large-Q 2 region, our calculations show that 
these form factors have about the same magnitude as the single-pion or two-pion analyses of the experimental data. 
However, in the low-Q 2 region, both the CLAS analysis and PDG value (at zero momentum transfer) still have 
opposite sign from ours for both proton and neutron transition form factors. 

In this workshop, we report preliminary results using Nf = 2 + 1 anisotropic dynamical lattice, focusing on the 
low-Q 2 region. Figure [2] shows three pion masses at a volume of 3, 2.5, 2.5 fm, where the volume dependence is not 
significant for Roper states; the experimental values are reconstructed using experimental helicity amplitudes from 
CLAS [8TfT2]. At the current statistics, the data indicate a sign change in the form factor relative to the quenched 
approximation, suggesting that in the low-Q 2 region, the pion cloud effects are significantly larger. Thus, removal of 
the quenched approximation causes a sign change. 

Axial Transition Couplings: The axial transition form factors contain spin-structure information about hadrons, 
and since the Roper is the lowest excited state of the nucleon, it is important to obtain this information for use in 
experimental analysis. Processes such as neutrino-nucleus scattering that are used to extract neutrino mass splittings 
and mixing angles require such inputs to control systematic errors. 

There are two approaches to extracting axial transition couplings on the lattice. A direct calculation of the transition 
form factor through the matrix elements 

(N^A^N^iq) = u N2 (p') 

at various momenta allows extraction of G\(q 2 ) and extrapolation to zero transfer momentum. We carry out the first 
lattice calculation of the axial proton-Roper transition form factor on the same 2+1-flavor anisotropic lattices at the 
heaviest pion mass of 875 MeV; see the left-hand side of Fig. [3j Our current estimation of G\ (0) is of the order of 0.1. 
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FIG. 3: (Left) Axial transition form factor G\ R (Q 2 ) on the N f = 2 + 1 anisotropic lattices with M w « 390, 450, 875 MeV. 
(Right) Pion form factor with multiple pion masses at 580, 875, 1350 MeV. The experimental points are shown as (black) circles 
while the lowest gray band is the extrapolation to the physical pion mass using our lattice points. 



Alternatively, we can study the volume-dependence of the nucleon mass, and use finite-volume chiral perturbation 
theory to extract the coupling constants. For example, NPLQCD studied the volume dependence of octet baryons 
at high precision. The axial decuplet-octet-meson couplings are extracted from finite-volume heavy-baryon chiral 
perturbation theory, |<7aat-7t | = 2.80(18)(21); similarly for other couplings such as |(7s* Att | = 2.49(23)(35), |3h*Htt| = 
2.49(23) (35). The advantage of the latter approach is that it is easier to calculate and obtains better signal-to-noise 
ratio than the direct approach. However, it depends on the reliability and our knowledge of the applicable chiral 
perturbation theory. A similar approach might be applied in the near future to find couplings such as Qa,n* an d 

9NN*tt- 

Large-Q 2 Prospects: Significant progress has been made in lattice QCD to achieve large-Q 2 form factors. Remain- 
ing obstacles are mainly due to the available lattice spacings 2 used in the current calculations and the signal-to-noise 
ratio at higher momenta for a given lattice spacing. Refs. E] investigate the latter limitation by reexamining the 
typical steps in a form-factor calculation: traditional methods optimize lattice parameters to eliminate excited-state 
contamination for a hadron at rest or low momentum. When one starts to boost the hadron to higher momentum, 
the hadron energy becomes compatible with the excited-state mass at rest, and the original parameters screen out 
any useful signal. The solution to this problem is to allow multiple parameter values that overlap with not only the 
ground-state hadron but also excited ones, and to extend the form-factor analysis to specifically include the excited 
states at each momentum. Thus, at higher momenta, there would be operators that still have good overlap. Using 
this approach, the number available form factors can be extended by a factor of 2-3, until the lattice discretization 
effects at 0((ap) 2 ) become dominant. The right-hand side of Fig. [3] displays an example using Nf = 2 + 1 pion form 
factors with 580, 875, 1350 MeV pion masses with Q 2 reaching nearly 7 GeV 2 [14 for the highest-mass ensemble. The 
extrapolated form factor at the physical pion mass shows reasonable agreement with JLab precision measurements. 
Future attempts will focus on decreasing the pion masses and exploring <5 2 -dependence of pion form factors for yet 
higher Q 2 . 

Refs. [HlHl] a l so suggest an approach to remove the first obstacle through a lattice technique called "step scaling", 
as demonstrated in Fig. 8 in Ref. 114} . using multiple lattice spacings and volumes to reach higher momenta by an order 
of magnitude. By calculating the step-scaling function at overlapping momentum points (or interpolating) we can 
reduce the systematic error due to finite- volume or lattice-discretization artifacts. However, generating several volumes 
and tuning for multiple lattice spacings of dynamical Nf = 2 + 1 lattices requires a large amount of computational 
resources; we hope such a proposal will become feasible as petascale computing facilities become available in the near 
future. 

Conclusion and Outlook: There has been much progress in Roper-related quantities calculated from lattice QCD, 
and the precision and techniques used have been greatly improved over the past couple years. We see impressive 
excited-baryon spectra from both the HSC and CSSM groups. The goals of the outlook from the previous Nstar 
Workshop, including the difficult particle identification from irreducible representations of the cubic group, have been 



2 The momenta allowed on the lattice (assuming periodic spatial boundary conditions) are ^o" 1 , where L is the length of the spatial 
volume, while n ranges 0, 1, v2, v3, etc. Thus, the finer the lattice spacing, the higher momentum transfer one can achieve. 
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implemented and published. Much remains to be done, but the progress already made has been remarkable. 

We see that finite-volume effects will play a significant role in the upcoming precision world. We note that even 
for the nucleon mass, the finite-volume effect is non-negligible for M V L w 4. With the Roper, the mass shift of the 
central values could be up to 170 MeV. We encourage the both HSC and CSSM to look for volume corrections in their 
calculations, either with multiple- volume calculations or estimating them from effective theory. A volume dependence 
study could also provide an indirect approach to extract the axial coupling constants related to Roper states. 

The Roper-nucleon electromagnetic and axial transition form factors also progressed to dynamical lattice calcula- 
tions. We use large volume anisotropic lattices to avoid finite-volume effects from the Roper state. The effort has 
been concentrated on the small-Q 2 region, where the sign discrepancy has been observed with respect to the previous 
quenched studies. We may have resolved the discrepancy in the small-Q 2 region from the CLAS analysis data that 
we have been seeing in quenched calculations. A first axial transition form factor is also reported in this proceeding, 
which will give us insight into spin structure. Further studies to increase statistics and refine the analysis on the 
current data are on-going. 
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